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Abstract 

A one-dimensional scattering problem off a 5-shaped potential is solved analytically and 
the time development of a wave packet is derived from the time-dependent Schrodinger 
equation. The exact and explicit expression of the scattered wave packet supplies us with 
interesting information about the "time delay" by potential scattering in the asymptotic 
region. It is demonstrated that a wave packet scattered by a spin-flipping potential can give 
us quite a different value for the delay times from that obtained without spin-degrees of 
freedom. 
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§1. Introduction and summary 

It is well known that the "time" occupies a rather peculiar position in quantum mechanics 
in the sense that it is just a parameter in the Schròdinger equation and no satisfactory "time 
operator" [1] has been found or devised so far, even though it is certainly measurable in 
experiments. This peculiarity can be one of the main sources of controvèrsies over, for 
example, the definition of the so-called "tunneling time" or the interaction time of a particle 
with a potential when the particle is represented by an almost monochromatic wave packet 
and its spatial width is longer than the range of the potential. 

To treat this kind of problem, one sometimes uses the time-índependent Schròdinger 
equation and forms a wave packet by superposing its solutions. Even though it can supply us 
with some information on the time development of the system under consideration, one has, 
at the same time, to remember in which sense this kind of treatment can represent the actual 
dynamical process: The solution to the time-independent Schròdinger equation is considered 
to describe the actual scattering process as an approximately stationary process, developing 
around the scatterer, where both the incoming and outgoing waves are present. The use of 
stationary solutions is justified in this restricted sense. Of course, we may construct a wave 
packet from these solutions with an appropriate weight function, in order to see the explicit 
time dependence. In this case, however, we have to be careful about the choice of the very 
moment t — since the actual scattering process is by no means translationally invariant 
(i.e. there is a moment when the particle is injected). In connection to this, we also need 
to correctly understand the meaning of the weight function. In this respect, it would be 
more natural and can be unambiguous to introduce a wave packet from the beginning to 
represent the incident particle and to examine its time evolution on the basis of the time- 
dependent Schròdinger equation. We shall set t = when the particle is injected and the 
weight function is so chosen as to reproduce the incident wave packet. 

The purpose of this paper is twofold. 

First, a one-dimensional time- dependent Schròdinger equation for a particle scattered by 
a simple 5-shaped potential is solved analytically to give an exact and explicit wave-packet 
solution in §2. Even though it remains at an elementary level, the derivation turns out to 
be rather instructive and suggestive. The exact expression of the wave packet enables us to 
understand how its asymptotic behavior emerges and is to be restricted on the energy shell 
as t — > oo. The connection to the stationary solutions is easily seen. The analysis is extended 
to incorporate possible internal degrees of freedom in §3. As a concrete and still solvable 
example, the scattering of a spin-1/2 particle off a spin-flipping 5-potential is considered and 
its wave-packet solution is derived. 

Second, in §4., we choose an almost monochromatic incident wave packet (i.e. spatially 
broad wave packet) and derive its asymptotic form, from which the mean position is calcu- 
lated and compared with that of a free wave packet without interaction. In the simplest case 
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with no internal degrees of freedom, it turns out that in spite of the vanishing interaction 
range of the 5-potential, the scattered wave packet exhibits a finite time delay or advance, 
depending on whether the interaction is repulsive or attractive. This delay (or advance) 
time is shown to agree with that derived from the energy derivative of the phase shift [2] 
caused by the potential scattering. In the spin-1/2 case, there are four scattering channels 
depending on whether the particle is transmitted or reflected and the spin is flipped or not. 
It is shown explicitly that contrary to the first case without spin-degrees of freedom, there 
remains no time delay or advance in any channels. This can be traced back to the reality 
of the scattering-matrix elements (i.e. the transmission and reflection coeficients) in this 
case. Notice that these two simple examples and their results imply that the notion of time 
delay of a wave packet can be quite different from the classical counterpart ~ [interaction 
range] /[particle's velocity] and its estimation fully needs quantum mechanical treatment. 

§2. Wave-packet solution for scattering off a 5-potential 

Let us consider a simple scattering problem, i.e. a one-dimensional scattering of a particle 
off a single 5-potential. The Hamiltonian is given by 

(2.D 

We shall solve the time- dependent Schròdinger equation for a wave function in momentum 
space (p\ij>)t, 

f) 2 

^<í#}í = (pm)t = W>* + gC(t), (2.2) 

where we have put 

The solution is easily found 

(p\^) t = e - r P 2t / 2mh (p\Íj) - i 9 - j* dt'e- ip2 ^l 2mh C{t') (2.4) 

in terms of the initial wave function {p\ip)o and C(t), the latter of which satisfies the consis- 
tency condition 

C (t) = í A e -^/2mft (g |^ )() _ £ í dq i* dt , e -i q \ t -t')/2mn c{tl y (2 5) 

J ItïTí n J 2nn Jo 

This equation can be reduced to an algebraic one by means of the Laplace transform 

C{s) = r dte' st C{t), (s>0). (2.6) 

J 



The solution C(s) is 



r (s) = 1 f d( i W)g (0 7) 

K) l + i{g/h)T{s) J 2Tïhs + iq 2 /2mh' K ' 



where ^ ^ 

In order to obtain C(t) from C'(s) via the inverse Laplace transformation, we need to make 
an analytic continuation of the latter from the original domain of s > into the left-half 
complex s-plane where Res < 0. Since the above function T is easily evaluated to be 



where y/s has a cut along the negative imaginary s-axis, we can write down C(t') (t' > 0) as 

ds 



dq , , (X f io °+ e ds Vïs{Vïs - i(g/Ti)^m/2Ti\ e st ' ( 2 ' 1 °' ) 
ioo+e 2tíí i(s — img 2 /2U 3 )(s + iq 2 /2mk)' 



Observe that the above integrand in (2.10) has a cut along the negative imaginary s-axis 
and two simple poles at s = —iq 2 /2mh (for fixed q) and s = img 2 /2h' i , the latter of which 
exists only on the second Riemannian sheet if the coupling constant g is positive (repulsive 
potential) and on the first sheet if g < (attractive potential). Therefore the original contour 
shall be deformed, avoiding these simple poles, to encircle in the left-half plane of the second 
Riemannian sheet. See Fig. I. 



Fig. 1. 



The choice of these deformed contours facilitates the integration over s and we obtain 

dq 



x 



\ e -iq 2 t'/2mh + \ e img 2 t'/2h 3 

l + img/k\q\ l + (hq/mg) 2 t 2 - 1 - 1 ) 

q írn^ f°° , ,/ye lyt ' 
+ -V< —V / dy- yy 



tiTi V 2h Jo (y - mg 2 l2T?) (y + q 2 /2mh) \ ' 

which is vàlid irrespectively of the sign of g. Each term in the square parentheses in (2.11) 
corresponds to the contribution arising from each singularity of the integrand mentioned 
above. 
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This explícit expression for C{t') leads us to an exact expression* of the wave function 
in p-space [see Eq. (2.4)] 



W)t = é 



-ip 2 t/2mh 



X 



+ img/h\q\ 

+ TTíW lAi[(p2 + (m9/fi)2)/2mR| 

+ JL flï-p f°° dy- ^ At ^ + p2 / 2mh ^ 



7rh V 2h Jo (y- mg 2 /2h 3 ) (y + q 2 /2m%) 



(2.12) 



where 



„ixt 1 

A t [x] = — . (2.13) 

ix 

Notice that this is an exact expression [3] for the wave function at t > and no assumption 
has been made in its derivation. 

Obviously the first term represents the free evolution of the incident wave packet, while 
the other terms stand for the scattered components. It would be interesting to observe 
that apart from the factor exp(—ip 2 t/2mh), the t dependence appears only through the 
functions A^s, each of which contains an energy-nonconserving component. The appearance 
of these energy-nonconserving terms is due to the energy-time uncertainty relation and they 
are expected to die out when the wave packet has passed through the interaction region 
and reaches the so-called wave zone. In order to understand this phenomenon clearly, let us 
consider the asymptotic limit t — > oo in (2.12). Notice that for large t, the function A t [x] 
approaches the Dirac 5-function 

A t [x] t ^2nS(x). (2.14) 

Therefore in this asymptotic limit, only one of three A^s, that is Af[(í> 2 — q 2 )/2mh] survives 
to give the energy conservation p 2 /2m = q 2 /2m, while the others A t [(p 2 + (mg/h) 2 )/2mh] 
and At[y + p 2 /2mh] vanish since their arguments are positive definite. It is worth noticing 
that the argument of At is nothing but the sum of the initial energy and the imaginary part 
of the location of the simple pole in the s-plane so that contributions coming from poles in 
the upper-half complex s-plane are to die out in the t — > oo limit, since they have positive 
imaginary parts and are not able to satisfy the energy conservation: Only those singularities 
located on the lower-half s-plane have relevance in the asymptotic limit. 



* It is possible to get an apparently different expression for {p\ip)t, if a different contour 
is chosen in (2.10), even though it is expressing completely the same quantity. The final 
form is dependent on the choice of the contour and one may choose a convenient contour, 
according to one's purpose. 
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It is now clear that in the asymptotic region t — > oo, the wave packet (p\ip)t is confined 
to the energy shell and is given by 



Observe that the two terms in the last expression correspond to the transmitted and reflected 
waves and the former is given by the sum of the free (nonscattered) and the scattered waves. 
Equation (2.15) shows clearly and explicitly the connection of the wave-packet (i.e. the 
time-dependent) solution to the stationary solution in the asymptotic region: In fact, the 
two factors in front of {p\ip)o and {—p\ip)o are nothing but the transmission (T) and reflection 
(TZ) coefficients, respectively, obtained for plane-wave scattering with defmite momentum p 

T(p) = - J—^ , Ttfp) = - Ímgl ^} í . (2.16) 
v ' l + img/h\p\ v ' l + img/h\p\ v ' 



§3. Scattering off a spin-flipping 5-potential and its wave-packet solution 

In this section, we consider a scattering of a particle endowed with internal degrees of 
freedom from a 5-shaped potential which is able to change the internal state of the particle. 
For definiteness and simplicity, let a spin-1/2 particle prepared in the up-state be injected 
at t — to the spin-flipping potential. Two spin-eigenstates of the particle, energetically 
separated by hu, shall be taken to be those of the third Pauli matrix 03, so that the total 
Hamiltonian is 

H=^- + ^(l + a 3 )+g5(x)a 1 . (3.1) 

The last term makes the particle's spin flip when it interacts with the potential. 

It is important to note that here we retain the spatial degrees of freedom of the particle. 
This is in contrast to the case in which, e.g. spin flip by a magnètic field is considered, 
where the spatial degrees of freedom are completely neglected and the interaction time is 
set to be [range of magnètic field] /[particle's velocity] from the outset. This is justifiable 
under the situation where the range of interaction is much larger than the width of the wave 
packet representing the incident particle, even if the latter is broad enough to approximate 
the particle as a plane wave. Obviously this condition does not hold in our case, since the 
interaction range of the 5-potential is 0. Actually, a question about the interaction time for 
such a short-range potential is one of the motivations of the present work. 

If we decompose the wave function in terms of momentum and spin eigenstates 

\f)t = Jdp [\p T) (p T \ri>)t + \P i) (P i W)t] , (3.2) 
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the Schròdinger equation becomes the coupled equations between (p j \ip)t and {p j \ip)t 



f) 2 

ih-^U^t^^ + ^ipU^t + gC^t), 



v 

2m 



(pimt + gC^t), 



where 



dq 



(3.3a) 
(3.3b) 

(3.4) 



Since the initial state is set in a spin-up state, (p [ |·0)o = and the solution to Eqs. (3.3) 
is given by 

(p | \j/;) t = e -i(p 2 /2ma+u;)í^ | _ ^ Ç dt < e -i(p 2 /2mh+u)(t-t') c ^ tl ^ ^ ^ 

n j o 

(p | = -i| J^dt'e-^^^C^t'). (3.5b) 

The consistency conditions that the above functions C|(í) and C|(í) should satisfy 

C Jt) = í ^ e -^ 2 /2rnn +í ,)t {q T |^ )o _ -| / ^ (ft / e -i(ç a /2 m a+ w )(t-Oc7,(í'), (3.6a) 
j 27ra « j 27m jo 



(3.6b) 



are converted to algebraic equations after the Laplace transformation. After an elementary 
calculation we find the solution 



cm 



-J 



dq 



(q T l^)o 



^s(s + iui) 



2nh s + i(q 2 /2mh + u) J s ( s + i w ) _ {i mg 2 /2T?) 



(g/k)Jim/2hs_ 



• (3-7) 



It is clear that the Laplace transforms Cj(s) and C^(s) have two simple poles at s = 
i{\joo 2 + (mg 2 /h 3 ) 2 — ui)/ 2 and s = —i(q 2 /2mh + lü), two branch points at s = and 
s = —íuj and a cut along the negative imaginary s-axis. The final expression of Cf(í') and 
Cj(í') is dependent on how the integration contour is deformed in the complex s-plane. 
For the evaluation of C|(í'), we shall deform the contour as in Fig. 2(a) 



Fig. 2. 
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to obtain 



x 



■ ™>g jtf /2mh+u)t' 

h 6 



^q 2 (q 2 + 2mhu) 



+ i 



g_ r^ v f 

kV 2h Jo 



dy 



^q 2 (q 2 + 2mhu) + (mg/K) 2 
Jyt' 



(3.8) 



7i y + {q 2 /2mh + u) Jy( y + üj- (mg 2 /2h 3 ) 



+ Íü?k) Jo 



dy 



-iyt' 



2ny- (q 2 /2mh + u) y(uj - y) + (mg 2 /2h 3 ) 2 

Inserting this into (3.5a) and integrating over t', we reach an exact formula for (p j \ip)t- 
In particular, in the asymptotic region t — > oo, the first, second and third terms in the 
above square parentheses give contributions proportional to 5(p 2 — q 2 ), S(y + p 2 /2mh + u) 
and S(y — p 2 /2mh — to), respectively, after the í'-integration, so that only the first term can 
survive to give the asymptotic behavior of (p f \ip)t 



(pUÏÏt t ^e- i <r 2 / 2mn+u '»(p J t\il>)o 

- 2 ' dq 



+ 



mg 
1? 



2nh 



q2 e -i(p 2 /2mh+uj)t 



I 

\lq 2 (q 2 + 2mhüü) + (mg/h) 



-A.T\mhb(p 2 — q 2 ) 



-i(p 2 /2mh+u)t 



v5V + 2mM 



\Jp 2 (p 2 + 2mhcü) + (mg/h) 2 
-(mg/h) 2 



(Pt 



+ 



p 2 (p 2 + 2mhuj) + (mg/hy 



{-p T 



(3.9) 



The first and second terms in the parentheses respectively stand for the components, trans- 
mitted and reflected through the potential without changing the spin state and the factors 
multiplying them coincide with the ordinary transmission (7+) and reflection (TZr) coefli- 
cients obtained by solving the time-independent Schròdinger equation. Observe that these 
coefficients are both real quantities. 

Similar treatment can be done for the evaluation of Cj(í'): We shall take the deformed 
contour, depicted in Fig. 2(b), and get 



c T (0 



dq 
2Ïrh 



(<zt 



q 2 (q 2 + 2mhu) 



q 2 (q 2 + 2mhüu) — (mg/h) 4 



,-i(q 2 /2mh+uj)t' 



— i- 



mg L 



h A 



v 



dy 



-iyt' 



>Jy(y-u) 



co 2rcy- (q 2 /2mh + cu) y(y - u ) - (mg 2 /2h 3 ) 2 _ 

(3.10) 



In the t — > oo limit, (p | |^)í is shown to behave like 



— — huj 

2m 



—ipt/2mh 



X 



p 2 (p 2 — 2mhüu) + (mg/h) 2 



(\/p 2 -2mhu t |^)o ( 311 ) 



Jp 2 (p 2 — 2mhuj) + (mg/k) 2 



These two terms correspond to the transmitted and reflected components of the particle, 
whose spin state is flipped to the down state through the interaction with the potential. The 
connection to the plane-wave solution is clear, since the common factor in the parentheses 
coincides with the transmission (TjJ and reflection (TZ\) coeficients, which are both real and 
take the same form in this particular case. Observe the presence of the #-function, which is 
necessary for energy conservation, since the total energy of the particle in the down state is 
given by p 2 /2m and the particle will gain energy heu by flipping its spin from up to down: 
No spin-down state with energy less than Heu is allowed because the incident particle has 
been prepared in the up state. 

§4. Estimation of the delay time 

The exact solutions to the time-dependent Schròdinger equations we have obtained in the 
previous sections explicitly show that their wave-packet solutions, even though they possess 
energy-nonconserving components within the finite time duration, are to be confined to the 
energy shell in the asymptotic (t — > oo) limit and converge to the states composed of the 
transmitted and reflected waves, multiplied by the relevant transmission (T) and reflection 
(1Z) coeflicients for the plane wave solution, that is, 

(p\i>) t *=3? e-tfW» [T(p)<pMo + K(p)(-p\m , (4.1) 

where E p is the energy of the state \p). For the moment, for notational simplicity, the 
dependence on the possible infernal degrees of freedom of the particle shall be suppressed. 

In order to illustrate the asymptotic behavior of the wave packet more clearly and quanti- 
tatively, let us express the above wave function in ordinary configuration space. We assume 
that the initial particle is represented by an almost monochromatic wave packet with a 
Gaussian profile 

( P - Pof 
A8p 2 h 1 

where M is a normalization constant and óp pq. This wave packet is distributed around 
its mean position xq < with width h/2Sp and moves with an average momentum po > 
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{p\^) =Me W \- KF - -(p- Po ) Xo \, (4.2) 



toward the potential. In order not to make any sensible overlap between the initial wave 
packet and the potential located at the origin, the inequality XQ + h/25p <C is also assumed. 

The transmitted wave packet in x-space is given by the Fourier transform of the first 
term of the RHS of Eq. (4.1) 



(x\ 



-fi=T(p)e- iE P t / n+tpx / n (p\^}o. 



(4.3) 



If the transmission coeficient T(p) is a slowly varying function of p, we may expand T(p) 
around po as a power series of p — po to perform a Gaussian integral over p in (4.3). Under 
this assumption, the transmitted wave packet (in the asymptotic region) turns out to take 
also the Gaussian form 



<a#)í r ~T(po) exp 
where 



| . . 2 % % -i 

— (x-x -pot/m-ihT'(po)/T(po)) + -p x - -E po t\, (4.4) 



Ah z u 



1 it 

U = T7. rT + 



45p 2 ' 2m%' (4 ' 5) 

It is straightforward, if lengthy, to extract the real part of the exponent in (4.4), which can 
be written in the following form 



5p 



x 



- x - vtrit - 5tr)] 2 + {àp 2 /pl) [ P oRe(T'( PQ )/T(p G 



h 2 i + ^E Sp t/hf 

where i>tr stands for the mean velocity of the transmitted wave packet 

1 + 2(5p 2 /p 2 )poRe(T'(p )/T(p ))' . 
The quantity <5t r in (4.6), given by 



(4.6) 



Vtr = — 

m 



(4.7) 



PO 



lm(T'( PQ )/T( P0 )) [l + 2(6p 2 /p 2 )p Re(r(p )/T(p c 



-i 



measures a temporal displacement of the transmitted wave packet and is called the time delay 
or advance, according to its sign. It is interesting to note that even though our derivation is 
quite different from that in [2], the main part of the above 5t T coincides with that obtained 
by taking the energy derivative of the scattering phase shift [2] : The correction, however, is 
due to the wave packet effect, which is obtainable only when the wave packet itself is treated 
directly. 

Similarly, we shall define the reflected wave packet (in the asymptotic region) by 



(x\ 



f ^L=TZ( P )e- iE ^ n+ipx / n (-p^) 



(4.9) 
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and perforin a Gaussian integral over p, under similar conditions as in (4.4). The behavior 
of the reflected wave packet, again shown to be a Gaussian, is read from the real part of the 
exponent 



5p 2 



n 2 1 + {AEspt/ny 



2 2 

x + xo + v ri (t - 5 I{ )] +(5p 2 /p 2 )[poRe{n'(p )/7Z( Po ))] , (4.10) 



from which we obtain its mean velocity and the time delay 



1 + 2(5p 2 /p 2 )p Re(n'(p )/n(p 



Po 

V T Ï = — 

m 

5 T{ = ^Im(^(po)/^(po)) [l + 2(5p 2 /p 2 )p Re(n\p )/n(p )) 



(4.11) 
(4.12) 



It would be interesting to estimate the value of the delay times for simple potential 
scatterings we have considered in the previous sections. Since the transmission and reflection 
coefficients for a particle with a definite momentum po, scattered by the 5-potential given in 
(2.1), are given by 

-iü(po) 



T(po) 



(4.13) 



1 + iíí(po) ' 1 + iíí(po) ' 

where íl(po) = m 9/hpo, the delay times for the transmitted and reflected wave packets are 



estimated to be 



5tr = ^9M l T {pQ )\ 2 \l + 2(5 P 2 / P 2 )\n(po)\ 2 ^ i 



2 E Po 
hQ(po) 



2 E, 



po 



\T( P0 )\ 2 l-2(5p 2 /p 2 ) \T( P0 



(4.14) 
(4.15) 



Notice that both delay times 5t T and S T { are positive, showing time delays of the wave packets, 
for a repulsive potential (g > 0), while, for an attractive one (g < 0), they become negative, 
which imply that the wave packets are actually advanced. Observe also that these time 
delays vanish not only for the weak coupling limit g — > 0, but also for the strong coupling 
limit g — > oo. Because in either limit, the wave packet is completely transmitted or reflected, 
we may understand that for the wave packet to exhibit a finite time delay or advance, the 
presence of both waves, transmitted and reflected, is crucial and that this phenomenon is 
ascribable to a kind of interference effect. This is still at a speculative level, however, we 
can say that the time-delay phenomena have an essentially quantum origin and can not be 
compared with the classical counterpart, which would also be clear from the fact that the 
5-potential has a vanishing range of interaction. 

The nontrivial nature of the time delay may be seen in the other case considered in §3., 
where a scattering of a spin-1/2 particle off the spin-flipping 5-potential has been analyzed. 
We have shown explicitly in (3.9) and (3.11) that the scattered wave packet in each channel 
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is represented just like (4.1) with channel-dependent transmission and reflection coeficients. 
If a particle prepared in the up state with definite momentum po is scattered by the spin- 
flipping potential given in (3.1), these coefficients read 



(4.16) 
(4.17) 




\JpHpI + Zmhuj) + (mg/h) 2 
.mg 1 



n ^JpUpI + 2mhw) + (mg/h) 




Notice that these quantities are essentially real and therefore no finite delay times are ex- 
pected for any channels. [See Eqs. (4.8) and (4.12).] This implies that no analogy to the 
previous simplest case without the spin-degrees of freedom is found! Even though the appear- 
ance or disappearance of such delay times in potential scatterings has yet to be understood 
well, we can conclude that its analysis certainly requires completely quantum mechanical 
treatment. 
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Figure captions 



Fig. 1. Singularities of the integrand in (2.10) and the deformed contours for (a) g > 
and (b) g < 0. There are two simple poles, denoted by crosses, at s+ = img 2 /2/ï 3 and 
s_ = —ip 2 /2mh and a cut, represented by a bold line, along the negative imaginary 
s-axis in either case. Dashed lines run on the second Riemannian sheet. 

Fig. 2. Deformed contours for the evaluation of (a) Cj(í') and (b) C-j·(í'). There are 
two simple poles, denoted by crosses, at s' + = t^uo 2 + (mg 2 /h 3 ) 2 - u)/2 and s'_ = 
—i(p 2 /2mh + a;), two branch points, denoted by dots, at s = and s = —íuj and a cut, 
represented by a bold line, along the negative imaginary s-axis. Dashed lines run on the 
second Riemannian sheet. 
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